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Abstract 

The moduli space of charge k SU (2) BPS monopoles is diffeomorphic to the moduli 
space of degree k rational maps between Riemann spheres. In this note we describe 
a numerical algorithm to compute the monopole fields and energy density from the 
rational map. The results for some symmetric examples are presented. 
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1 Introduction 



The BPS monopole fields, $, Ai, are SM(2)-valued fields in satisfying the Bogomolny 
equation 

Di(^ = --EijkFjk (1.1) 

and the boundary condition that = — |tr$^ = 1 at infinity. The Higgs field, $, 

defines a map from the two-sphere at infinity to the two-sphere of vacua and the winding 
number, k, of this map is a non-negative integer which equals the magnetic charge of the 
monopole (in suitable units). 

The moduli space of charge k monopoles, Aik, is a 4/c-dimensional manifold and it has 
been known for some time |Q that Aik is diffeomorphic to the space of degree k based 
rational maps between Riemann spheres. More recently, a new correspondence between 
Aik and the space of degree k unbased rational maps of the Riemann sphere has been 
proved. In this new approach the rational map arises as the scattering data, along half- 
lines from the origin, of a linear operator constructed from the monopole fields. In this 
note we describe a numerical inverse scattering algorithm to compute the monopole fields 
and energy density from the rational map. The results of the algorithm are displayed for 
several symmetric examples. 

The algorithm involves solving a nonlinear elliptic partial differential equation which is 
equivalent to the Bogomolny equation, but for which the boundary conditions are given in 
terms of the rational map. It is the ability to specify a unique solution through the rational 
map boundary conditions which is the crucial feature here. It eliminates a problem which 
would arise in attempting a numerical solution of the original Bogomolny equation, in that 
it is not at all clear how to specify the boundary conditions on the monopole fields so as 
to select a unique monopole from the 4/c-dimensional family of solutions. 

An alternative numerical method for constructing monopoles, based on the Nahm trans- 
form, already exists and it is perhaps useful to make a brief comparison of the two 
approaches. The numerical Nahm transform requires the solution of only ordinary dif- 
ferential equations and not a partial differential equation, so from that point of view it 
requires less computational resources. However, the input required is Nahm data, which 
involves the analytic solution of a matrix ordinary differential equation. The general so- 
lution, for charge k, involves theta functions on a curve of genus {k — 1)^ and this is too 
difficult to obtain explicitly in practice. For certain symmetric examples p|, |^ the curve 
has symmetries which allow the solution to be obtained on a quotient curve which is ellip- 
tic, and it is such examples which can then be used as input for the numerical code. Some 
special non-elliptic examples have been studied numerically by using approximate Nahm 
data 1 12], but even these cases rely on very symmetric configurations and there is currently 
no method available to deal with the general case for charge greater than two. In contrast 
the input for the numerical algorithm presented in this note is a rational map and so the 
data is free. Any rational map can be chosen as the input data and hence the code can be 
used to obtain the monopole fields for all points in and not merely some special cases. 

Let us briefly review the main points of p| to see how the rational map arises as scatter- 
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ing data. First introduce polar coordinates, r, 6, (p, and combine the angular coordinates 
into the Riemann sphere parameter z = e*"^ tan(^/2). The rational map is obtained by 
considering Hitchin's equation 

{Dr-i^)s = (1.2) 

for the complex vector s, along each radial half-line from the origin out to infinity, with 
the direction of the half-line determined by the value of z. The boundary condition on 
$ implies that there is a one-parameter family of solutions s{r) = {wi{r) , W2{r)y that 
are bounded as r oo. Define R to be the ratio of these components at the origin ie. 
R = Wi (0)7^2(0). It can then be shown [^] that i? is a holomorphic function of z with 
degree equal to the monopole charge k. Thus R : (DP^ 1— > (DP^ is the required degree k 
rational map between Riemann spheres. Note that the rational map is unbased, since 
a gauge transformation replaces R by an SU{2) Mobius transformation determined by 
the gauge transformation evaluated at the origin. Thus the correspondence is between a 
monopole and an equivalence class of rational maps, where two maps are equivalent if they 
can be mapped into each other by a reorientation of the target Riemann sphere. Since the 
above construction does not break the 5*0(3) rotational symmetry of P^ then a monopole 
which is invariant under a subgroup G C 50(3) will have an associated rational map R 
which is G-invariant (up to Mobius transformations). 

2 The Algorithm 

Writing the Bogomolny equation ( |1 . 1| ) in terms of the coordinates r, z, z and choosing 
the (complex) gauge 

$ = -iAr = --R-^drH, A, = R-^d^H, = (2.1) 
where H G SL{2, (D) is a Hermitian matrix, results in the single equation |p 

dr (h-^ drH) + ^^^jf^' d, [H"' d,H) = 0. (2.2) 

In P it is proved that solutions of this equation are uniquely determined by boundary 
conditions at large r which are in one-to-one correspondence with rational maps between 
Riemann spheres. This result is presented more explicitly in where it is shown that 
the connection between the rational map R{z) and the matrix H{r, z) is the large r 
asymptotic relation 

The input to the algorithm is the rational map i?, and the main task is to compute the 
solution of equation (|2.2|) subject to the boundary condition that as r ^ 00 the solution 
has the asymptotic form (|2.3|). 
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To compute solutions of the elliptic equation ( p.2| ) we use a standard heat flow approach 
and study the parabolic equation 



H-^ dtH = dr {h-' drH) + ii±m_5^_ (^-1 d^H) . (2.4) 
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As t ^ oo the solutions of this equation are static, that is t-independent, and hence solve 
the original elliptic equation. 

Although the choice of coordinates, r, z, z, has played a crucial role in the analysis so 
far, these are not good coordinates in which to implement a numerical solution. In R'^ the 
most efficient coordinates to use for a numerical code are Cartesian coordinates, in which 
equation (|2^ ) becomes 



Ht = Hu - ^XiHi - HiH-'Hi + -SijuXiH^R-^ H^. (2.5) 



Here an index i denotes partial differentiation with respect to Xj, and the summation 
convention is used over the spatial indices. 

Although we are now using Cartesian coordinates the method by which we have ob- 
tained equation (|2.5|) is reflected in the broken translational symmetry. Equation ( p.5|) can 
not be applied at r = 0, but this is not a problem since there is an additional constraint 

that at the origin H = I, the 2x2 identity matrix. 



We solve equation ( p.5| ) using a flnite difference approximation, with spatial derivatives 
approximated by second order symmetric differences and the time derivative approximated 
by a flrst order forward difference. The discretization is performed on a regular grid, 
typically using 70^ grid points. On the boundary of the grid the matrix H is flxed to be 
the asymptotic expression ( ^73|) . Note that it is a subtle matter to enforce the boundary 
condition numerically since the elements of the matrix H are unbounded as r oo. Thus 
the size of the bounding box must be big enough to ensure that the monopole conflguration 
sits comfortably inside the box, but the bounding box must not be too big otherwise the 
unbounded growth of the elements of H will lead to numerical difficulties. In practice we 
find that a bounding box with |x| < 5 is satisfactory for all the examples we present. 

The initial conditions are taken to be 

f 2g{r) (\R\^-l -2R \\ 
^ = '"nTT^l -2R l-\R\^)j (2.6) 

where g{r) is a monotonically increasing function with g{0) = and g{r) ~ r for large r. 
Note that from the relation (|2.1| ) between H and $ this initial condition gives a Higgs fleld 
whose length, is spherically symmetric. There are no spherically symmetric SU{2) 

monopoles with > 1, so it is clear that the initial condition is not a good approximation 
to the flnal solution for charges greater than one. 

The elements of H are evolved directly in the time evolution and so the unit determinant 
constraint is preserved by performing the rescaling transformation 

H , , 

E ^ , 2.7 
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after each iteration. A typical example requires around 5000 iterations to converge to a 
static solution. 

From the solution H the monopole fields are computed from the relations ( |2.1j ), where 
derivatives are again approximated by second order symmetric finite differences. For ex- 
ample, the length of the Higgs field, which is a gauge invariant quantity, is given by 

||$||2 = --tT{drHdrH-^). (2.8) 

8 

The energy density, S, is then obtained from the Higgs field via the formula 

S = dA\\n'. (2.9) 

In Figure 1 we display the results of our algorithm, in the form of energy density isosur- 
face plots, for several examples. We have chosen some particularly symmetric examples, 
with charges from two to seven, and the explicit rational maps for each example are listed 
in Table 1. 



Charge 


Rational Map 


Symmetry 


Figure 


2 




Axial 


la 


3 


z''' — \Piiz 


Tetrahedral 


lb 


4 


2*-2V3i22+l 


Cubic 


Ic 


5 


2^—52 

1-52-* 


Octahedral 


Id 


6 


4 1 ■ 

2 +ja 

ja2'^+2^ 


Dihedral 


le 


7 


2^-725-722-1 
2^ + 725-722 + 1 


Dodecahedral 


If 



Table 1 : Rational maps for the monopoles shown in Figure 1. 

The details of these rational maps and a discussion of their symmetries can be found in 
ref.P]. A reason for choosing the examples with Platonic symmetry is that these cases fall 
into the special class which have been computed with the numerical Nahm transform 0, ^ 
and hence we can verify that the results obtained from the two very different numerical 
approaches are in agreement. The charge six example has only a dihedral symmetry, 1^4^, 
and the associated Nahm data is not known, since it is not an elliptic example. This case 
is therefore one in which the numerical Nahm transform could not be used to construct the 
monopole and hence this is the first computation of this particular monopole. The degree 
six rational map given in Table 1 has Dn^d symmetry for all values of the real parameter 
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a, and the particular example shown here is for a = 0.16. The reason for this particular 
choice is that there is an ansatz to obtain approximate Skyrmions from rational maps 
0, and within this ansatz the rational map with a = 0.16 is the one which minimizes the 
energy of the approximate Skyrmion. Of course, for monopoles all configurations of a given 
charge have the same energy, but it is interesting to note that the monopole energy density 
isosurface in Figure le is remarkably similar to the Skyrmion baryon density isosurface 
shown in ref. It has already been observed that the known Platonic monopoles resemble 
the corresponding Skyrmions but a deep understanding of this fact is not yet available. 
Within the rational map approximation for Skyrmions there is an understanding of the 
relationship between the Skyrmion baryon density and the rational map and it would be 
interesting if a similar relationship could be found to connect the monopole energy density 
to the rational map. 



3 Outlook 

We have presented an algorithm which computes monopoles from rational maps and 
displayed the results for some examples. In this section we comment on some possible 
applications of this code. 

One of the surprising results obtained from the numerical computation of monopoles is 
the discovery that a charge k monopole can have more than k zeros of the Higgs field |^ ^ . 
The total number of zeros counted with multiplicity is k, but some of the zeros can have 
a negative multiplicity, which are known as anti-zeros. The current numerical evidence is 
obtained using the numerical Nahm transform and it is satisfying to confirm the existence 
of anti-zeros with the algorithm presented in this note, since the two numerical approaches 
are very different. Anti-zeros are indeed detected in the new algorithm, as demonstrated 
in Figure 2. In Figure 2 we plot the length of the Higgs field, along the xs-axis 

for the octahedral 5-monopole shown in Figure Id. The solid curve displays the initial 
condition, where the only Higgs zeros is at the origin, whereas the dashed curve represents 
the solution at the end of the iteration. It can be seen that there are now three zeros along 
this line, and hence by the octahedral symmetry there are seven in total. The zero at the 
origin is an anti-zero giving a total multiplicity of five as required. At present there is 
little understanding of the anti-zeros phenomenon and no signature for their appearance 
is known in terms of the associated rational map. However, now that an algorithm is 
available to construct a monopole from its rational map it will be possible to make a more 
extensive study and investigate any possible connections with the topology of the space of 
rational maps. 

In testing conjectures of electric-magnetic duality it has proved extremely useful to use 
the correspondence between monopoles and rational maps |11|. Some preliminary attempts 



0] have been made to use the rational map correspondence to determine which monopole 
configurations play a key role in this context and perhaps the construction of the monopole 
fields themselves may prove useful. 

In the geodesic approximation [|TU| the dynamics of k slowly moving monopoles can be 
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approximated by geodesic motion in the fc-monopole moduli space. It is possible to find 
geodesies as fixed point sets of group actions on the space of rational maps. By constructing 
the monopole fields for such families of rational maps it would be possible to study various 
monopole scatterings in some detail. 

Finally, the method introduced here could be extended to construct SU {N) monopoles 
and this may be of use in understanding non-abelian clouds 0. 
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Figure Captions 

Figure 1: Energy density isosurfaces for the monopoles listed in Table 1. 

Figure 2: Plots of ||$||^ along the xs-axis for the octahedral 5-monopole. Initial condi- 
tion (solid curve) and final solution (dashed curve). 
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This figure "monratfigl.gif" is available in "gif" format from: 



http://arXiv.org/ps/hep-th/9905066vl 



